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1. INTRODUCTION 

Today, Doubly Fed Induction Generator (DFIG) is widely used for power generation in wind farms 
because of its high efficiency, energy quality and the ability to control power supplied to the network [1-3]. 
However, DFIG-based wind turbines have a number of disadvantages such as their sensitivity to load 
disturbances, turbine’s rotational speed variation, and internal and external parameters of the system. 

Similarly, the system is multi-variable, dynamic, highly coupled and non-linear, making control 
very difficult [4]. On the other hand, the need for development of wind farms for the use of renewable 
sources, pushed the evolution of research in the field to improve the efficiency of the electromechanical 
conversion and the quality of the electrical energy supplied by different controls. [5]. The main idea of this 
approach is to analyze the stability of the nonlinear system without solving the differential equations of this 
system. It is a very powerful tool to test and find sufficient conditions for the stability of different dynamic 
systems. The study presented in this paper aims to optimize the energy performance of a wind turbine in 
order to maximize the wind energy captured while reducing structural problems cited previously and improve 
control performance by reducing time of calculation. This optimization is capable to minimize the costs of 
electricity production. 
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In this paper, we will start by exposing a wind turbine modeling. Then, a study of the operating 
point tracking technique and maximum power will be presented. Subsequently, we present a DFIG model in 
the "dq" frame and the general principle of controlling both power converters which is based on the control 
strategy (FOC, and Sliding Mode). 

In this work, we develop two control strategies entitled: Field Oriented Control (FOC), and Sliding 
Mode Control. A study of these commands is detailed and developed in Matlab / Simulink environement. 
Finally, the results are analyzed to compare the effectiveness of the controlled system of the two 
control strategy. 


2. WIND TURBINE MODELLING 
2.1. Energy conversion chain 

The overall diagram of a wind energy conversion chain connected to the power grid is depicted 
in Figure 1. 
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Figure 1. Diagram of an energy conversion chain 


2.2. Modelling of the wind turbine Power 
In this part, we define the mathematical model used in modeling the wind system. The expression of 
the mechanical power captured by the wind turbine and transmitted to the rotor (Rotor) is expressed by [3]: 


Protor = Pwina- Cy (A, B) 


1 
Pwind = > PSV* 


Protor = > PS. Cp(A, B)v? (1) 


With Cp(A, B)is the power coefficient, v is the wind speed (m/s), p is the air density and S is the 
swept area (m?) 

Betz proved that the maximum extractable power for a given wind speed by ideal turbine and 
conditions is Cpmax (à B) = 0.593, this limit is known as the Betz limit. For the DFIG, it is possible to 
model the power coefficient (C, (A, 6)) with a single equation that depends on the tip speed ratio (A) and the 
pitch angle (B) of the blades. This can be written as: 


1 T1 

C(A, f) = cy. (c2.2— c3.B - ca). FA + Cé.À (2) 
With 
= 1.5872 
c, = 116 
C3 = 0.4 
G= 
Cs = 21 
c, = 0.0085 
1 1 0.035 
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j= (3) 


V 


The wind turbine is characterized by the curves of the coefficient power, which corresponds to the 
aerodynamic torque (T;-o¢9,), Shown in following equation. 


Protor 1 1 
Trotor — a = 5 pS. C(A, pv? (4) 


It is necessary to connect the wind turbine and the DFIG to the Gearbox in order to adapt the speed 
of the turbine to that of the generator. It is modeled by the two following equations: 


Q=G. Q 
Trotor = G. Tm (5) 


The evolution of the generator speed (Q) depends on the fundamental equation of the dynamic that 


characterizes the mechanical behavior of the turbine and generator from the total mechanical torque applied 
to the rotor. It is given by the following formula (Figure 2): 


aQ 
JT = Tn — Tem — f.Q (6) 





L'arbre 





Initialise wares Turbine 
Figure 2. Wind-turbine model 


2.3. Mechanical regulation of the power of a wind turbine-pitch control 

A wind turbine is sized to develop on its shaft a power called Pn nominal power obtained from a 
wind speed vn, called nominal speed. When the wind speed is greater than vn, the wind turbine must modify 
its parameters in order to avoid mechanical destruction, so that its rotation speed remains practically 
constant. [6-10]. 

Beside the nominal speed v,,, we also specify: 
a. The starting speed vq, from which the wind turbine starts to supply energy 
b. The maximum wind speed v,,, for which the turbine no longer converts wind energy, for reasons of 

operational safety. 

The speeds v, ‚Vq and Vm define four zones on the diagram of the useful power as a function of the 
wind speed: 
Zone I, where P = 0, the turbine does not work. 
Zone II, in which the power supplied on the tree depends on the wind speed v. 
Zone III, where the speed of rotation is kept constant and the power P supplied remains equal to Pn. 
Zone IV, in which the operating safety system stops the transfer of energy [6-10]. 

So, the control of the wedge angle is in zone III in order to maintain, using the orientation of the 

blades of the turbine, the power extracted at a value called nominal value 

We present below the influence of the variation of the angle of calage on the value of the power 
coefficient, which represents the energy efficiency of the turbine (Figure 3 and Figure 4) 


aoo op 
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Figure 3. Diagramme de la puissance utile sur l’arbre Figure 4. Influence of the angle B on the 
en fonction de la vitesse du vent power coefficient 


It can be seen that as the wedge angle increases, the efficiency of the turbine decreases. 
The control of the wedge angle $ for the orientation of the blades of the turbine can be achieved at a 
single loop control or two (Figure 5 and 6): 





Figure.6. Control of the angle B using double-loop control 


The control performed by a single loop control is used to control the nominal power. The output of 
the regulator of the nominal power then corresponds to the reference value of the angle B, and it is 
transmitted to the system of orientation of the blades by an open-loop control. 

When the control is carried out with double-loop control, the first corresponds to the control of the 
nominal power, and the second to the control of the calibration angle. 

The second control mode is more complex to implement in terms of the synthesis of regulators, and 
it is more expensive since it requires the use of two regulators for its control. 


2.4. Principle of power control 

There are two principles of aerodynamic control to limit the power extracted from the turbine to the 
rated power of the generator: [2-6] 

a. "pitch" or "variable pitch" system that adjusts the lift of the blades to the wind speed to maintain a 
substantially constant power in the speed zone II 

b. "stall" or "aerodynamic stall" system, the most robust because it is the shape of the blades that leads to a 
loss of lift beyond a certain wind speed, but the power curve drops faster: it is therefore a 
passive solution. 

The first principle is associated a blade orientation mechanism, allowing the variation of the 
wedging angle during the operation of the wind turbine to allow it to adapt to different wind conditions. The 
interest of this control appears by observing the characteristics of Figure 7 that present the power of the 
turbine as a function of the speed of rotation for different wind speeds. 


Int J Pow Elec & Dri Syst Vol. 10, No. 4, Dec 2019 : 2101 — 2117 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 2105 


Puissance mechanique [Watts] 
Vitesse cu vent [m/s] 





Vitesse de rotation de la génératrice | tr/mn] 


Figure 7. Power of a turbine according to its speed of rotation, Set in wind speed. 


2.5. MPPT strategy 

We have seen above that the extraction of wind power depends on the value of the power 
coefficient. The figure, below, represents the power according to the rotational speed of a given 
turbine (Figure 8): 

It can be seen that, for each value of wind speed V, the characteristics power - rotational speed 
graph goes through a maximum corresponding to a maximum power reached by an optimal rotational speed. 
However, the maximum power is obtained for a maximum power coefficient (Figure 9): 
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Figure 8. MPPT-Power of wind turbine Figure 9. MPPT-Power coefficient cp 


This maximum value of the power coefficient corresponds to an optimal specific speed / opt. This 
specific speed is obtained for an optimal rotational speed: 


Acc 
ont = a (7) 





This optimal value of the rotational speed is sent as a reference for the speed control of 
electrical machines. 

The speed control of the electrical machines can be achieved either in scalar control or in vector 
control (also called Field Oriented Control). The scalar command is used for the steady-state control while 
the vector control is used in transient mode. As part of this thesis, we opted for vector control, since the 
models used are developed under transient conditions. The control is more precise and faster, and moreover it 
allows a control of the magnitudes in amplitude and phase. 

Also, the control of the static converter in the machine side makes possible to control the torque of 
the generator in order to obtain the desired rotational speed to extract the maximum power. This consists in 
controlling the electrical machine by a general cascade structure using nested loops. 

Torque and flux, therefore, will be controlled by very fast internal current loops. These commands 
are carried out in a rotating frame (vector control with flux oriented). The torque setpoint is obtained from a 
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slower external speed control loop. As for the speed reference, it is calculated from a MPPT strategy as 
presented above. 


3. MODELLING OF THE DFIG 

The generator chosen for the conversion of wind kinetic energy is the double-feed asynchronous 
generator [11-12]. 

The modeling of the DFIG is described in the Park reference system. The following equation system 
describes the global modeling of the generator that connects the voltages to the currents, the speed, the flow 
of the machine, and also the mechanical equation. 

3.1. Electrical equations 

Stator and rotor voltage 


apy 
Vas = Rslas + a = Os Pqs 
dp 
Vas = Rel gs + ae — O;Pas 
apa 
Var sat Ry lar + EA HE 0 Par 
dPar 
Vor = Ry lor ar a oi 0 Par (8) 


Stator and rotor currents 


i 1 M 
ds o. L Pas D. L Par 
M 

= Pas TP 
2 Wde S Ddue 
Pae 1 M 
dr — D. iz Par D. L Pas 

1 M 

lar = Fi Var x Pas (9) 
Stators and rotor flux 
Pas — Lilas + Ma, 
Pas = Lslgs + MI qr 
Par = Lyla, + Mas 
Par = Lrlgr + Mis 
Tem =T, + fX+J— (10) 


Stator and rotor powers 

The control of the DFIG must allow an independent control of the active and reactive powers by the 
rotor voltages generated by an inverter. 

In the two-phase reference system, the active and reactive stator and rotor powers of an 
asynchronous generator are written: 


3 

R= z (Vas- las + Vs) 
3 

Qs = z (Vas- las + Vale) 


3 
P, = 2 (Var- lar + Velay) 
3 
Q; = = (Var: ar + Viele) (11) 


3.2. Mechanical equations 
The electromechanical torque can be determined by the following equations (Figure 10): 
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Figure 10. DFIG model 


Tem = PUas: Pas — lgs: Pas) (11) 


For obvious reasons of simplification, a d-q reference linked to the stator rotating field and a stator 
flux aligned with the d axis have been adopted [11]. In addition, the stator resistance can be neglected since 
this is a realistic assumption for generators used in wind power. On the basis of these considerations, the 
torque becomes 


M 
Tem = P T, (Pas: lar Fe Paslar) (12) 


4, FIELD ORIENTED POWER CONTROL 
4.1. Principe of power control 

The principle of stator flux orientation control is based on [13]: @g, = Ps and Pgs = 0. 

For large wind turbines, the stator resistors are negligible, and also the flux becomes constant. For 
this, the flux has an expression: 


Pas = OS Ll gs + M. or 
Pas — Lslas + M.lar (13) 


From stator and rotor equations, we obtain: 


M 
las = Le 
1 
las = 1, (Ps z M.lar) (14) 


The expression of the electromagnetic torque [14]: 
M 
Tem = PPs: las = —P 7, Oster (15) 


Based on the assumptions: 


Vas = 0 
Vas me = Ws. Ps 
Ps = p, (16) 


The active power P, and the reactive power Q, of the wind system can be written [14] 
Fe — Vas dás + Vas: las 





Q; = Vas- las = Vas: ls (17) 
From the previous equation, we obtain [14]: 
E; = Vs. los = Vs- 7- lar 
S 
Vs? M 
Qs = Vo-las = E- Vee lar (18) 
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The rotor flux is given by: 





=|L u lar + M Ms 
Par = r L; "idr wsLs 


Par — (L, — or (19) 


Subsequently, the control vectors Vg, and V,, can be determined by the equations: 


M? M? 
Var — Q +S. (L; = D) «dar — Ws. g. (z = X) Lor 
S S 


M? M? M 
Var = |(Ry +S. (Ly | Igy — @s-9- (Ly - ©) Tar + 9. Ve (20) 


The following figure shows the diagram block using Simulink that models the equations 
defined previously: 

The Figure 6 shows the modeling of DFIG equations using first-order transfer functions for 
voltages. The vector control will be easily implemented later, considering that the slip value g is negligible 
and the influence of the coupling is weak. 


4.2. FOC technique 


In this work, we use two types of vector control. The first is based on controlling the rotor voltage of 
the DFIG to control the power of the system (Figure 11). 








PID Cont oir 











Figure 11. Direct vector control 


The second type (indirect control) is based on controlling the rotor current to control the power. The 
control vector is generated from Equations 19 and 21 (Figure 12). 


4.3. Simulation results 


In this part, we present the simulation results of the proposed model for direct and indirect control of 
flux (FOC) (Figure 12). 
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Figure 12: FOC control 








Figure 13. Direct vector control of a wind system based on DFIG 


The Figure 14 present the simulation results for the proposed control models: 
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Figure 14. Results of simulation of a wind system based on the DFIG variable speed 


5. SLIDING MODE POWER CONTROL 
5.1 Sliding mode power control strategy 


The principle of this technique consists in bringing the states trajectory of a system towards the 
sliding surface and to switchs it by means of appropriate switching logic around it to thesequilibrium point 


(Figure 15); hence, the phenomenon of slip [15-16]. 





Figure 15. The trajectory in the phase plane 
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The design of the sliding mode control algorithm is carried out mainly in three complementary steps 
defined by [15-16]: 


5.1.1. Sliding surfaces: 
Considering the following non-linear system: 


x = f(x,t) + B(x,t).u(x, t) (21 


Where: 

xER” : Statut Vector, f (x, t)eR” , B(x,t)eR™™” , ueR™ : Statut Vector 

In order to ensure the convergence of a state variable to its reference value, it is necessary to choose 
a sliding surface which is a scalar function,so that ,the variable to be adjusted slides on the latter. 

J. SLOTINE proposes a form of general equation [15-17] given by: 


S(x) = (= + 5)" .e(x) (22) 


Where: 

6: Positive gain ; (x) = x, — x: Deviation of the variable to be regulated; n: is a relative degree 
which represents the number of times that the surface must be derived to make the command appear . 

The objective of the command is to keep the surface at zero. The latter is a linear differential 
equation whose unique solution is e(x) = 0. 


5.1.2. Conditions of existence and convergence: 

In order to ensure the convergence mode, LYAPUNOV proposes a function V (x) which guarantees 
the stability of the nonlinear system and the attraction of the variable to be controlled to its reference value, 
the latter is defined as follows: 


V(x) ==.S(x)? > 0 (23) 


For decreasing the function (x) , it suffices to ensure that its derivative V(x) = S(x).S(x) < 0. The 
idea is to choose a scalar function S(x) to guarantee the attraction of the variable to be controlled to its 
reference value, and to design a command" U '",, so that the square of the surface corresponds to a function of 
LYAPUNOV [16]. 


5.1.3. Determination of the control law 
Once the sliding surface and the criterion of convergence have been chosen, it remains to determine 
the control necessary to attract the state trajectory towards the surface and then towards its equilibrium point 

while maintaining the conditions of existence of the sliding mode. 
The command (uw) is a variable structure command given by [18] 

+ 

= a s(x,t) >0 (24) 
u(x), s(x,t) <0 


Where: u ‘and u~are continuous functions, with w= # ut 

This control (wu) of discontinuous nature will force the trajectories of the system to reach the sliding 
surface and to remain there in the vicinity of the latter despite the presence of the disturbances. 

The system trajectories on the surface S are not defined since the command (u) is not defined for 
S = 0, to do this, FILLIPOV [19] and UTKIN [20] propose a method called « The Equivalent Command » 
[16, 21] (Figure 16). The addition of this command pre-positions the system in a desired stable state [22] and 
in addition reduces the CHATTERING phenomenon. 


Robust power control methods for wind turbines using DFIG-generator (Imane El Karaoui) 


2112 O ISSN: 2088-8694 


Perturbation 


SONY Switching Law 





Figure 16. The equivalent command 


This structure consists of two parts [16], one concerning the exact linearization (Ueq) and the other 
the stability(U,). 


U = Usa + Un (25) 


(Ueq is used to keep the variable to be controlled on the sliding surface (x) = 0 . It is deduced, 


considering that the derivative of the surface is zero S (x) = 0. 
Considering the state system (17), we sought to determine the analog expression of the control U. 
The derivative of the surface S (x) is given by the following equation: 


as as as 

S(x) = 5, £0) +5 BW). Ueq + 5. BX). Un (26) 
In the sliding and permanent regime, the sliding surface is zero and its derivative as well as the 

discontinuous part are also null, so the expression of the equivalent command becomes as follows: 


ane FOG oe) (27) 


. ð . ox 
It is necessary that = .B(x) #0 , so that the equivalent command can take a finite value. 


By replacing U,,with its expression inS (x), we obtain: 
as 
S(x) = 5, BUX): Un (28) 


In order to satisfy the condition of attractiveness S(x). S (x) < 0, the sign U,,of must be opposed to 
that of =. B(x). S(x). 

The discontinuous control (U,,) forces the dynamics to converge towards the surface and ensures the 
insensitivity of the system with respect to uncertainties and perturbations. 

Several forms are proposed in the literature [21, 23], the simplest one is given by the following 
expression U,, = K. sign(S(x)) 

Where: K is a positive constant and (sign) the classical sign function. 

But the main disadvantage of the relay type control is the phenomenon well known by 
"CHATTERING". In steady state, the latter appears as a high frequency oscillation around the equilibrium 
point (Figure 17), because of the very discontinuous nature of the sign function. 

To remedy this problem, several investigations have been made. One of the solutions envisaged 
consists in introducing a stop band around the switching surface. To do this, it suffices to replace the function 
(sign) with the saturation function (Sat), whose discontinuities in the vicinity of zero are less brutal. 

This saturation function can be expressed by the following equation [15-16]. 


lifo>e 
sat(~) -1 if p <E (29) 
zif lel <e 
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S(x)=0 
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Path / 


Figure 17. Phenomenon of "CHATTERING" 


5.1.4. Application of the sliding mode command to the DFIG. 

After presenting the theory of Sliding Mode Control, we will analyze in this part the application of 
the sliding mode control to the DFIG. 

According to the Field Oriented Control principle, the quadrature component of the stator field is 
forced to zero and the direct component is equal to the total stator field. 

The preceding equations become as follows: 


























Va = Ù 
Vsq = Vs = Ws. Os 
vra = Rp. lyg + ri + Wy. Brg 
Vq = Rr- lrg +E + op. Dra (30) 
pP e= “Veoh 
s 

payta 

SL SL T4 
Ing = = 7 = | rq a Sma (31) 


5.2 Active power control surface 
The sliding surface proposed by JSSLOTINE is given by the following equation: 


S(x) = (=< + 5)" e(x) (32) 


To control the active power, we take n = 1, the expression of the active power control surface is 
as follows: 


S(P;) = e(P;) = Psref — Es (33) 
Its derivative is given by equation (32): 
S(P,) = Piref -P (34) 


We replaceP, by its expression as well as Lg obtaining us: 





i ‘ M Vr Ry MVs 
SPs) = Psref + V-T- C m a. — Wr. ) (35) 


Lao. Lio. 4 Ls.Lr.0.Ws 


By replacing the expression of with , the command clearly appears in the following equation: 





5 . M Vr e +Vz n Rr MVs 
SE) = Péref + Ve — ee 6 ) (36) 


be Ly Leo Ebel pO We 
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In sliding and permanent mode: 





S(P)=0 

S(P)=0 

Van = 0 

Vigeq = —Ls bp. a Pore + (Ree Log OME) (37) 


The discontinuous control is given Vg, by: 
Van = Ks. sat (S(P;) (38) 
Where:K,: Positive gain 


5.3 Reactive power control surface 
We take the same surface as that of the active power: 


S(Qs) = e(Qs) = Osref — Qs (39) 


The derivative of the surface is: 
S(Qs) = Or =- Q; (40) 


We replaceQ,by its expression as well asl,,, we get: 








S(Qs) = Qsres + Ven GE. La) (41) 


Ly. Ly.o | 


By replacing the expression of V,, with V,,4eq + V,un we obtain the following equation: 





; . M /VrdeqtVran Rr 
S(Qs) = Qsref + Vea a = Laird — Wr. Irq (42) 
In sliding and permanent mode: 
5(Q,) = 0 
S(Qs) = 0 
Vran =0 (43) 


Hence, the formula of V,ge,becomes as follow: 


Vrdeq = —Ls L L 


Lee Qsref + Rr- lra — Le. O- Trg. Wr (44) 


The discontinuous control V,g,1s given by: 


Vran = Ka- sat(S (Qs) (45) 
Where: K,: Positive gain 


5.4 Simulation results: 

In order to validate the robustness of the sliding control applied to a wind system based on the 
DFIG, the following figures show the results obtained from this application, where the following 
observations can be distinguished: 
— The wind profile is given by (Figure 18). 
— The aerodynamic power according to the MPPT has the same shape as that of the wind profile. 
— The power coefficient Cp and the specific speed are kept around their optimal values 0.5505 and 8 

respectively, which ensures maximum mechanical power 
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The rotation speed of the machineis less oscillating with a high reliability with respect to the 


input speed. 


The shape of the electromagnetic torque of the DFIGis negative with a shape different from that of the 


wind profile this is due to the dynamic torque of the inertia. 


Figure 18 shows the active stator power and its reference profile injected into the network, where a very 


good decoupling can be seen between the active and reactive power of the stator. 


The stator reactive power and its reference profile are shown in Figure 18, which gives a unit power 


factor factor cos @ (@) = 1 because (= 0) .Consequently, iron losses are minimized. 


Figure 19 shows the stator currents that are purely sinusoidal with a constant frequency 50 Hz and 


consequently a minimization of harmonics and therefore a better injection into the grid. 
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Figure19. Indirect field oriented control Result 
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In this section ,an analysis and interpretation of the results of the application of the sliding control to 
a wind system based on the DFIG will be presented as well as a comparison with a literature based on the 
same control [24] . 

The figures obtained present the different performances of the wind system, which show the 
robustness and good tracking of the references. The specific speed and the power coefficient vary in 
proportion to the wind speed, which shows the reliability of the system. The active and reactive power 
follows their references well, making the system more robust. 

In paper [24], it is clear that the active and reactive power exhibit many oscillations, which implies 
less robustness. It can also be noted that the frequency of the stator current is different from that of the grid, 
which is one of the major problems of wind systems. By comparing these results with ours, we note that the 
wind system in our work is more reliable and robust, which makes it possible to make a better injection of 
power in the grid. 


6. COMPARISON BETWEEN DIFFERENT CONTROL STRATEGY: FOC AND SLIDING 
MODE 
6.1. Field oriented power control 
Advantage: 
— The switching frequency is fixed; 
— The harmonic content of currents is well defined; 
— There is an application of the eight voltage vectors that can supply the voltage inverter; 
— Good quality of regulation of the current in steady state; 
Disadvantages: 
— The general structure of the control algorithm is complex to implement; 
— Transient dynamics are slow; 
— The parameters of the control algorithm depend on the parameters of the DFIG. 


6.2. Sliding mode power control 
Advantage: 
— The general structure of the control algorithm is simple to implement; 
— Robustness ; 
— The parameters of the control algorithm are independent of the parameters of the machine; 
— Very good dynamics during transient regime; 
Disadvantages: 
— The zero voltage vectors are not applied. 


7. Conclusion 

This work was devoted to modeling, simulation and analysis of a wind turbine operating at variable 
speed. The use of the two controls linear and non-linear (FOC, Sliding Mode) with the application of the 
MPPT and Pitch Control techniques, allows to improve the system performance and guarantee a good 
monitoring of the active and reactive powers. 

The comparative study, between these two controls, shows that the Sliding Mode control offers 
several advantages such as the robustness and the linearization of the system than the other commands. 
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